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Abstract 

How to identify the influential spreaders in social networks is crucial for acceler¬ 
ating/hindering information diffusion, increasing product exposure, controlling dis¬ 
eases and rumors, and so on. In this paper, by viewing the k-shell value of each node 
as its mass and the shortest path distance between two nodes as their distance, then 
inspired by the idea of the gravity formula, we propose a gravity centrality index 
to identify the influential spreaders in complex networks. The comparison between 
the gravity centrality index and some well-known centralities, such as degree cen¬ 
trality, betweenness centrality, closeness centrality, and k-shell centrality, and so 
forth, indicates that our method can effectively identify the influential spreaders in 
real networks as well as synthetic networks. We also use the classical Susceptible- 
Infected-Recovered (SIR) epidemic model to verify the good performance of our 
method. 
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developing efficient strategies to control epidemic spreading, accelerate infor¬ 
mation diffusion, promote new products, and so on. In view of this, many 
centrality indices have been proposed to address this problem, including de¬ 
gree centrality [12], betweenness centrality [13], neighborhood centrality [H] 
and closeness centrality [15], etc. In particular, Kitsak et al. proposed a /c-shell 
decomposition method to identify the most influential spreaders based on the 
assumption that nodes in the same shell have similar influence and nodes in 
higher shells are likely to infect more nodes, k-shell method is found to be 
better than the degree centrality index in many real networks [T]. However, 
recent researches have demonstrated that the nodes within the same shell often 
have distinct influences, and this method may fail in some networks without 
core-like structure, e.g., Barasasi-Albert network [16]. Thus, after this, some 
methods were proposed to further improve the performance of the k-shell 
method. For example, Zeng et al. proposed a mixed degree decomposition 
method by incorporating the residual degree and the exhausted degree HHi 
Liu et al have demonstrated that the existence of the core-like groups can 
result in the invalidation of k-shell method [18], and then they showed that 
the accuracy of k-shell method can be improved once the redundant links in 
networks are removed [13]. Chen et al. designed a semi-local index by con¬ 
sidering the next nearest neighborhood [2D]; Lin et al. presented an improved 
ranking method by taking into account the shortest path distance between 
a target node and the node set with the highest k-core value [21]; Recently, 
Bae et al. dehned a novel measure-coreness centrality index, which is given 
by summing all neighbors’ k-shell values [22] . 


In general, a node’s influence is not only dependent on its nearest neighbors 
but also on the nodes who are not the nearest neighbors [23E3] , meanwhile, 
their interaction influence commonly decreases with their shortest path dis¬ 
tance. If the k-shell value of each node is viewed as its mass, and the shortest 
path distance between two nodes is dehned as their distance, then we can use 
the idea of gravity formula proposed by Isaac Newton to measure the inhuence 
of nodes. Inspired by these factors, in the work, we propose a new centrality 
index to measure the inhuence of nodes, which is called gravity centrality in¬ 
dex. We apply the susceptible-infectious-recovered (SIR) spreading dynamics 
to evaluate the ehectiveness of our proposed method, the experimental results 
indicate that gravity centrality index can better evaluate the inhuence of nodes 
than the ones generated by degree centrality, betweenness centrality, k-shell 
centrality, closeness centrality, and so on. 


The layout of the paper is as follows: In Sec. [21 we hrst briehy review several 
typical centrality indices which are used to compare in this work, and the 
description of our method is presented. Then the experimental results are 
presented in Sec. [3l Finally, Conclusions and discussions are summarized in 
Sec. m 
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2 Method 


An undirected network is represented by G = {N, M) with N nodes and M 
edges, and its structure can be described by an adjacent matrix A = 
where = 1 if node i is connected to node j, and aij = 0 otherwise. 

Here we briefly review the definitions of several centrality indices that will be 
discussed in this work. 

The degree centrality (DC) of a node is defined as the number of nearest 
neighbors. The betweenness centrality (BC) of a node is defined as the frac¬ 
tion of all shortest paths travel through the node. The closeness centrality 
(CC) of a node is defined as the reciprocal of the sum of the lengths of the 
geodesic distance to every other node. The k-shell decomposition method (ks) 
is implemented by the following steps: Firstly, remove all nodes with degree 
one, and keep deleting the existing nodes until all nodes’ degrees are larger 
than one. All of these removed nodes are assigned 1-shell. Then recursively 
remove the nodes with degree no larger than two (i.e., remove all nodes with 
degree two, and keep deleting the existing nodes until all nodes’ degrees are 
larger than two.) and include them to 2-shell. This procedure continues until 
all nodes have been assigned to one of the shells HT]. 

To improve the exactness of k-shell method, the mixed degree decomposition 
(MDD) method was proposed by Zeng et al. [1^. The mixed degree km{i) for 
a node i is dehned by considering the residual degree kr{i) and the exhausted 
degree ke{i) simultaneously, which is written as: 

kraii) = Kii) + X*ke{i). (1) 


At each step of the MDD procedure, the nodes are removed according to the 
mixed degree, and the mixed degrees of remaining nodes are also updated. 
Where A is a tunable parameter between 0 and 1. As in Ref. nn, we take 
A = 0.7 in this work. 

Recently, Bans et al. designed a ranking method-neighborhood coreness Cnc 
by considering the degree and the coreness of a node simultaneously, the Cnc{i) 
for a node i is defined as [22] 

C„,(i) = E (2) 

jeA, 


where Aj is the neighbor node set of node i. They further developed an ex- 
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tended neighborhood coreness C'„c+, which is described as: 


= Y. CncU). 

jeAi 


( 3 ) 


Chen et al. proposed a semi-local centrality measure as a tradeoff between 
low-relevant degree centrality and other time-consuming measures (labeled as 
SL index). It considers both the nearest and the next nearest neighbors. The 
semi-local centrality SL{i) of node i is dehned as [20] 


Q(i!) = i: N{j), 

(4) 

j&As 


sL(i) = o(s), 

(5) 


sGAi 


where A* is the neighbor node set of node i. N{j) is the number of the nearest 
and the next nearest neighbors of node j. 


It is fact that, on one hand, the influence of a node is increased if its neighbors 
(here the neighbors of a node do not just includes its nearest neighbors, which 
may also include next nearest neighbors, next-next nearest neighbors, etc.) 
have higher value of ks] on the other hand, the interaction effect between two 
nodes decreases with their distance. Enlighten by the idea of classical gravity 
formula proposed by Isaac Newton, we can view the k-shell value of node i 
as its mass, and the shortest path distance between two nodes in network is 
viewed as their distance. In this way, the influence of node i is measured by 
(labeled as G): 


G{t) 


E 


ks{i)ks{i) 

4i 


( 6 ) 


where dij is the shortest path distance between node i and node j. is the 
neighborhood set whose distance to node i is less than or equal to a given value 
r. To reduce the algorithm complexity, in the paper, we let r = 3, i.e, only 
nearest neighbors, next nearest neighbors and the next-next nearest neighbors 
are considered. 

An extended gravity index is further developed based on Eq. ([6]), which is 
dehned as (labeled as G+): 

G+{t) = GU), ( 7 ) 


Aj is the nearest neighborhood of node i. 


4 



3 Experimental results 


In this section, we compare the effectiveness of other indices with onr G or G+ 
index from different aspects and on different networks, inclnding real networks 
as well as synthetic networks. 

We employ the standard susceptible-infected-removed (SIR) model [25] to es¬ 
timate the real spreading influence of the nodes (labeled by R). In detail, to 
check the spreading influence of one given node, we set this node as an in¬ 
fected node and the other nodes are susceptible nodes. At each time step, each 
infected node can infect its susceptible neighbors with infection probability (3, 
and then it recovered from the diseases with probability fi. In this paper, we 
set /i = 1.0. This process repeats until there has no any infected nodes. At last, 
the number of recovered nodes is used to reflect the real influence of the node. 
To guarantee the reliability of the results, all of them are at least averaged 
over 1000 independent realizations. 

First, a small example network with A^ = 20 nodes and M = 25 edges is given 
in Fig. 1 to intuitively compare these indices, the ranking lists from different 
indices are presented in Table 1. From Table 1, one can see that the ks central¬ 
ity cannot well distinguish the influence of nodes. Even in the same shell, the 
nodes’ influence may be totally different. Moreover, the result indicates that, 
our proposed G and G+ index can effectively identify the influence of nodes, 
i.e., the ranking list determined from G or G+ index is in good agreement with 
the ranking list obtained from SIR model in the last row. 

Table 1 

The ranking lists determined by different indices. Degree centrality: DC; mixed 
Degree decomposition: MDD; gravity centrality: G; extended gravity centrality: G+; 
extended neighborhood coreness defined in Eq. ([3|): C„c+; k-shell decomposition: ks; 
betweeness centrality: BC; closeness centrality: CC; semi-local centrality measure: 
SL; the node spreading influence evaluated by SIR model: R, by taking j3 = 0.25. 


Rank 

DC 

MDD 

G 

G+ 


ks 

BC 

CC 

SL 

R 

1 

1, 2 

1, 2 

2 

2 

0 

0, 2, 3, 4 

2 

0 

0 

2 

2 

0,4 

0, 4 

0 

0 

2 

1, 8, 10, 12, 14-16 

0 

2 

2 

0 

3 

3, 8, 10, 14 

3 

4 

4 

4 

others 

1 

4 

4 

4 

4 

6, 12, 15, 16 

8 

3 

3 

3 

— 

4 

1 

3 

1 

5 

others 

10,14 

1 

1 

1 

— 

14 

3 

1 

3 

6 

— 

12, 15, 16 

8 

8 

8 

— 

6, 8, 10 

8 

8 

8 

7 

— 

6 

14 

— 

10 

— 

16 

14 

10 

14 

8 

— 

others 

10 

— 

12, 14 

— 

15 

10 

— 

10 

9 

— 

— 

— 

— 

— 

— 

others 

— 

— 

— 


To validate the effectiveness of the G or G+ index, we apply it to 9 real net¬ 
works, including Facebook (Slavo Zitnik’s friendship network in Facebook) [26] . 
Netsci (collaboration network of network scientists) [23) Email (e-mail net- 
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Fig. 1. An example network with 20 nodes and 25 edges. Nodes with larger degrees 
have larger size. 

work of University at Rovira i Virgili, URV) [2S], TAP(yeast protein-protein 
binding network generated by tandem affinity purification experiments) [T7] . 
Y2H (yeast protein-protein binding network generated using yeast two hy¬ 
bridization) [29], Blogs (the communication relationships between owners of 
blogs on the MSN (Windows Live) Spaces website) [30], Router (the router- 
level topology of the Internet) [31], HEP (collaboration network of high-energy 
physicists) [32], PGP (an encrypted communication network) [33]. For simplic¬ 
ity, these networks are treated as undirected and unweighed networks in this 
work. The detailed information about these 9 real networks are presented in 
Table 2. 

How to improve the resolution is the key issue of an algorithm, for instance, 
in Ref. [M], Zhou et al. have clarihed that the resolution problem is a major 
reason for the poor performance of common neighbor index subject to the 
AUC value in link predication, and then they proposed a local path index to 
solve the resolution problem. Similarly, a good index in ranking the influences 
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Table 2 

Basic structural properties. N and M are the number of nodes and edges, respec¬ 
tively. f3th is the epidemic threshold. H is degree heterogeneity, given by {k‘^)/{k)‘^. 
f is assortativity coefficient. C is clustering coefficient. L is average shortest path 
length. D is diameter. 


Network 

N 

M 

(3th 

H 

r 

G 

L 

D 

Facebook 

324 

2218 

0.047 

1.567 

0.247 

0.465 

3.054 

7 

Netsci 

379 

914 

0.125 

1.663 

-0.082 

0.741 

6.042 

17 

Email 

1133 

5451 

0.053 

1.942 

0.078 

0.220 

3.606 

8 

TAP 

1373 

6833 

0.061 

1.644 

0.579 

0.529 

5.224 

12 

Y2H 

1458 

1948 

0.140 

2.667 

-0.209 

0.071 

6.812 

19 

Blogs 

3982 

6803 

0.072 

4.038 

-0.133 

0.284 

6.252 

8 

Router 

5022 

6258 

0.072 

5.503 

-0.138 

0.012 

6.449 

15 

HEP 

5835 

13815 

0.110 

1.926 

0.185 

0.506 

7.026 

19 

PGP 

10680 

24316 

0.053 

4.147 

0.238 

0.266 

7.463 

24 


of nodes should also has a high resolution. As illustrated in Table 1, G or G+ 
index is good at distinguishing the nodes’ difference, which is much better 
than the ks index. So to quantitatively measure resolution of different indices, 
a monotonicity index M{X) for a ranking list X is used [22]: 


M(X) 


M 1) -|2 

^ N{N - 1) ^ 


( 8 ) 


where N is the size of network, and Nf. is the number of nodes with the same 
index value c. If M{X) = 1, which means that the ranking method is perfectly 
monotonic and each node is categorized a different index value; otherwise, all 
nodes are in the same rank as M{X) = 0. The monotonicity M for different 
ranking methods is summarized in Table 3. Generally, the results suggest that 
G or G+ index can give higher value of M. Moreover, M{G) and M(G+) are 
very near 1 in some networks. Therefore, gravity method can better distinguish 
the node’s influence than other indices. 

The Kendall’s tan rank correlation coefficient r is used to measure the cor¬ 
relation one topology-based ranking list and the real spreading capability R. 
Let (xj, Ui) and {xj,yj) be a randomly selected pair of joint observations from 
ranking lists X and Y, respectively. If one has Xi > Xj and yt > yj or Xi < Xj 
and yi < yj, the observations {xi,yi) and {xj,yj) are said to be concordant. If 
Xi > Xj and y^ < yj or Xi < Xj and y^ > yj, they are said to be discordant. If 
Xi = Xj or yi = yj, the pair is neither concordant nor discordant [3331] • t is 
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Table 3 


M (.) is the monotonicity of the corresponding measures. 


Network 

M(DC) 

M(MDD) 

M(G) 

M(G+) 

M(Cnc+) 

M(ks) 

M(BG) 

M(GC) 

M(SL) 

Facebook 

0.9315 

0.9729 

0.9999 

0.9995 

0.9995 

0.8445 

0.9855 

0.9953 

0.9999 

Netsci 

0.7642 

0.8215 

0.9949 

0.9951 

0.9893 

0.6421 

0.3387 

0.9928 

0.9939 

Email 

0.8874 

0.9229 

0.9999 

0.9999 

0.9991 

0.8088 

0.9400 

0.9988 

0.9999 

TAP 

0.8991 

0.9599 

0.9994 

0.9994 

0.9981 

0.8380 

0.9238 

0.9988 

0.9992 

Y2H 

0.4884 

0.5304 

0.9966 

0.9960 

0.9633 

0.2972 

0.5063 

0.9957 

0.9936 

Blogs 

0.5654 

0.5906 

0.9976 

0.9976 

0.9868 

0.4670 

0.4004 

0.9973 

0.9971 

Router 

0.2886 

0.3009 

0.9967 

0.9965 

0.9657 

0.0691 

0.2983 

0.9961 

0.9953 

HEP 

0.7654 

0.8314 

0.9998 

0.9999 

0.9917 

0.6303 

0.5651 

0.9998 

0.9990 

PGP 

0.6193 

0.6678 

0.9995 

0.9997 

0.9851 

0.4806 

0.5099 

0.9996 

0.9986 


defined as 


N^-N2 
0.5N{N -1)' 


where A^i and N 2 are the nnmber of concordant pairs and discordant pairs, 
respectively. 

When we employ SIR model to check the spreading inflnence of nodes, the 
infection probability f3 shonld not be too small or too large. The epidemic 
cannot successfnlly spread over networks if {3 is too small, so the spreading 
capability of each node cannot be measured. On the contrary, if (3 is too large, 
the epidemic can easily outbreak over almost whole network, leading to the 
spreading capability of each node cannot be distinguished too. Thus, in this 
work, we hrst obtain the epidemic threshold fith for each network, which is 
given as Pth ~ (^)/(^^)) with {k) and be the average degree and the 
second order average degree ESI, respectively. The value of Pth for different 
networks is given in Table 2 too. Then, we choose the value of P to be slightly 
larger than the threshold Pth when computing r for different indices (a new 
index to measure the influence of nodes was proposed in Ref. [36], which is 
independent on the parameter P). The results in Table 4 manifest that our 
method outperforms the other methods in most cases. 

To further estimate how the infection probability P affects the effectiveness 
of different methods, the correlation value r as a function of P for different 
methods is shown in Fig. 2. As described in Fig. 2, in most cases, G or G+ index 
provides better performance than the other index when P > Pth (the values of 
Pth for different networks are illustrated by the dot lines in Fig. 2). However, 
Fig. 2 clearly indicates that though the global indices, such as betweenness 
index and closeness index are time-consuming, they are not good at measuring 
the influence of nodes in these networks. Meanwhile, the performance of MDD 
method in identifying the node’s influence is almost the same as the degree 
centrality. 























Table 4 

r(.) is correlation of corresponding methods for given /3. 


Network 

/3 

tdc 

'I'MDD 

'^G 

TG+ 


'^ks 

tbc 

'J'CC 

'^SL 

Facebook 

0.050 

0.767 

0.796 

0.861 

0.913 

0.916 

0.735 

0.364 

0.720 

0.940 

Netsci 

0.130 

0.599 

0.620 

0.830 

0.852 

0.847 

0.525 

0.308 

0.330 

0.806 

Email 

0.070 

0.771 

0.790 

0.887 

0.937 

0.935 

0.779 

0.625 

0.822 

0.935 

TAP 

0.065 

0.725 

0.746 

0.870 

0.899 

0.873 

0.690 

0.273 

0.527 

0.886 

Y2H 

0.160 

0.445 

0.463 

0.827 

0.833 

0.825 

0.407 

0.412 

0.701 

0.775 

Blogs 

0.075 

0.525 

0.532 

0.834 

0.763 

0.795 

0.482 

0.390 

0.579 

0.706 

Router 

0.075 

0.322 

0.323 

0.797 

0.805 

0.786 

0.186 

0.315 

0.642 

0.790 

HEP 

0.110 

0.487 

0.506 

0.787 

0.865 

0.735 

0.485 

0.345 

0.784 

0.840 

PGP 

0.055 

0.479 

0.490 

0.784 

0.770 

0.756 

0.439 

0.313 

0.636 

0.747 


Previously, our results were obtained by setting the value of r = 3 (i.e., only 
the effects of the nearest neighbors, next nearest neighbors and the next-next 
nearest neighbors are considered). To check the sensitivity of r on our results, 
the effect of the parameter r on the value of r is plotted in Fig. 3. As shown 
in Fig. 3, in generally, the optimal value of r is about 3-5, and the value of r 
becomes stable when r is further increased. As a result, it is unnecessary to 
choose too large value of r, which just increases the algorithm complexity of 
our method. 

Besides the real networks, we also compare the performance of our method 
with other methods on the two typical synthetic networks-Barabas-Albert 
(BA) networks [16] and the Watts-Strogatz (WS) small-world networks [37] 
with N = 1000. Starting from a connected network with mo nodes to con¬ 
struct a BA network, at each step, a new node is added to the network and 
connected to m existing nodes according to the preferential attachment mech¬ 
anism, where m < mg [IE]. We set the number of nodes mg = 10 in this 
paper. The WS small-world model considers a ring nearest neighbor coupled 
network with N nodes. Each node symmetrically connects to its 2K nearest 
neighbors. Starting from it, a fraction p of edges in the network are rewired, 
by visiting all K clock-wise edges of each node and reconnecting them, with 
probability p, to a randomly chosen node [37]. During the rewiring process, 
self-connection and reconnection are forbidden. 

For BA network (see Fig. 4 (a) and (b)), one can see that the performances 
of G, G+ and Cnc+ indices are almost the same. The reason is that the three 
indices are all the improved methods of k-shell method, however, all nodes 
in BA network are almost classihed into the same shell when using the the 
k-shell method (so we do not calculate the case of ks in Fig. 4). Moreover, the 
results show that the three indices are better than CC index and are much 
better than DC, BC and MDD indices. For WS network (see Fig. 4 (c) and 
(d)), whose degree distribution shows Poisson distribution, i.e, their degrees 
are not so different. In this case, it is difficult for DC index to distinguish the 
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Fig. 2. (Color online) The value of r obtained by comparing the ranking list gen¬ 
erated by the SIR model and the ranking lists generated by the topology-based 
method on Facebook, Netsci, Email, TAP, Y2H, Blogs, Router and HEP. The dot 
lines correspond to the epidemic threshold. 

influence of nodes. However, as shown in Fig. 4(c) and (d), as /3 > f3th, the 
performances of G and G+ indices are still better than the other indices. In 
particular, for WS network, one can observe that the performances of G and 
G+ indices are much better than the Gnc+ index when {3 > (3th- The results in 
Fig. 4 suggest that our method can not only identify the influential nodes on 
real networks but also on synthetic networks. 


4 Conclusions and discussions 


In summary, in this paper, we have proposed a gravity method to identify 
the influential spreaders in complex networks. In the model, each node’s k- 
shell value is considered as its mass and the shortest path distance between 
two nodes is viewed as their distance. The idea of the gravity method comes 
from the well-known gravity formula, which is very dramatic and impressive. 
What’s more, the gravity model can reflect the facts that, on one hand, the 
interaction influence between two nodes is proportional to their corresponding 
k-shell values; on the other hand, the influences of the neighbors decreases with 
their distance. We employed our method on some real networks and synthetic 
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Facebook 
•— Netsci 



2 4 6 8 10 12 14 16 18 20 22 24 26 

r 


Fig. 3. (Color online) The effect of considered range r on the Kendall’s tan rank 
correlation coefficient r. Here the value of /3 for each network is the same to the value 
of fi in Table 4. We should address that, though the average distance of Facebook 
is about 3, the value of r is larger than 3, since the distance between a pair of 
randomly chosen nodes may be larger than 3. Therefore, the largest value of r is 
the diameter of the network. For example, r = 7 for the Facebook network. 

networks, by calculating the monotonicity index M, we found that our method 
can better distinguish the difference of node influence than other indices. Also, 
by computing Kendalls tan rank correlation coefficient r, we have shown that, 
in most cases, our method has a better performance in evaluating the node’s 
influence than other indices. Therefore, our method provides an effective way 
to identify the influential spreaders in social networks. 

Some extensions may be made based on this method. For example, by dehning 
the combination of node’s degree and node’s strength as the weighted degree 
of a node in weighted networks, Garas et al have proposed a new k-shell 
decomposition method for weighted networks [38]. Therefore, once the new 
k-shell value for each node in weighted network is assigned, our method can 
be simply generalized to weighted networks [39]. Also, if we view the closeness 
centrality, degree centrality, eigenvector centrality, and so forth as the mass of 
a node, then the gravity method may be further generalized. 

We only investigated the performance of the gravity method in some typical 
networks, and the classical SIR model was used to mimic the the spreading 
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p 



Fig. 4. (Color online) (a) BA: m = 3; (b) BA: m = 4;(c) WS: K = ?>, p = 0.05; (d) 
WS: K = 4, p = 0.05. The pink dot lines correspond to the epidemic threshold. 

dynamics. In reality, the strnctnre of networks and the spreading dynamics are 
diverse. For example, recent researches have illustrated that networks in na¬ 
ture do not act in isolation, but instead exchange information and depend on 
one another to function properly, that is to say, natural systems are organized 
in interconnected networks ^011411142] ; And some real spreading dynamics like 
the diffusion of rumors or opinions, the rise of scientihc ideas |^44II45] are 
different from the spreading of epidemic, which may challenge the effective¬ 
ness of proposed indices. For example, Borge-Holthoefer et al. have stated 
that the influential nodes in networks are absent when considering the rumor 
dynamics [1], also in Ref [5], authors have illustrated that the roles of nodes 
are dependent on the collective dynamics. Therefore, feasible methods need 
to be examined, we here hope our work inspire possible solutions to the above 
mentioned problems in the near future. 
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